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ABSTRACT 

We consider the curvature emission properties from relativistic particles streaming 
along magnetic field lines and co-rotating with pulsar magnetosphere. The co-rotation 
affects the trajectories of the particles and hence the emission properties, especially the 
polarization. We consider the modification of the particle velocity and acceleration due 
to the co-rotation. Curvature radiation from a single particle is calculated using the 
approximation of a circular path to the particle trajectory. Curvature radiation from 
particles at a given height actually contains the contributions from particles streaming 
along all the nearby field lines around the tangential point, forming the emission cone of 
1/7. The polarization patterns from the emission cone are distorted by the additional 
rotation, more serious for emission from a larger height. Net circular polarization can 
be generated by the density gradient in the emission cone. For three typical density 
models in the form of core, cone and patches, we calculate the polarization profiles 
for emission generated at a given height. We find that the circular polarization could 
have a single sign or sign reversal, depending on the density gradient along the rotation 
phase. The polarization profiles of the total curvature radiation from the whole open 
field line region, calculated by adding the emission from all possible heights, are similar 
to that from a dominating emission height. The circular polarization of curvature 
radiation has sign reversals in the patchy emission, while it has a single sign for the 
core emission, and is negligible for the cone emission. 
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1 INTRODUCTION 

Pulsar radio emission is generally believed to be generated 
by relativistic particles streaming out along the open mag- 
netic field lines in pulsar magnetosphere. The observed radio 
emission are generally highly li nearly polarized, and have 
significant circu l ar polarization (JLvne fc Manchester! 119881 : 
iHan et al.l 119981 : iRankin fc Ramachandranl 120031 ). Various 
radio emission mechanis ms are suggested to explain po- 
larization profiles (e.g. Be skin et al.l llQSSi : IXu et al.l I2OO0I : 
lGangadharall2010l ). 

Curvature radiation serves as one of the most pos- 
sible emission mech a nisms for pulsar radio emission. 
iBuschauer fc BenfordI (|l976l ) firstly developed the general 
formalism for coherent curvature radiation by a relativistic 
plasma streaming along curved trajectories using antenna 
mech anism, assuming the p l asma is perturbed by a plane 
wave. iBenford fc Buschauerl (|l977l ) applied this mechanism 
for pulsar radio emission and calculated the luminosity. 
lOchelkov fc Usovl (1980) analytically calculated the pulsar 
emission spectrum from curvature radiation. In polarization 
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aspect, it is known that c urvature radiation is alm ost com- 
pletely linearly polarized. iGil fc Snakowskil l|l990r) pointed 
out the sign rev ersal feature f or cir cular polarization (CP) 
of core emission. ICangadharal l|2010l ) considered the geome- 
try of the emission region in pulsar magnetosphere, obtained 
the polarization states self-consistently, and explained the 
correlation bet ween the positio n angle (PA) and the sign 
reversal of CP pan et al.lli998l ). 

Because pulsar magnetosphere is co-rotating with the 
neutron star, the streaming particles should have co-rotation 
velocity additional to the streaming velocity along the mag- 
netic field line, which affects t he curvature radiation prop- 
erties. iBlaskiewicz et al.l l|l991f) first proposed a relativistic 
model for pulsar polarization by incorporating the rotation 
effect, and predicted the phase lag between the centers of the 
PA curve and the profile. The rotation also influences the 
inte nsities of the leading and trailing components of the pro - 
filc ('Blaskicwic z et al.lll99il : iThomas fc Gangadharall2007l ). 
Dyks ct al. (2010|) analyzed the rotation and found two com- 
peting effects: the enhancement of the trailing side caused 
by aberration and retardation (AR) and the weakening of 
the trailing side caused by the rotation-in duced asymmetry 
of curvature radius. In the pioneer paper. IBlaskiewicz et al.l 
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l|l99ll ) calculated only the emission of particles with veloc- 
ity pointing towards the observer, and did not consider the 
emissions from particles of nearby field lines, which may still 
beam towards us. Meanwhile, the influence of ro tation on po- 
laxization intensities are not taken into account. iGangadharal 
l|201Cl ) considered the emission from all nearby field lines, 
but did not incorporate the rotation in his detailed work. 

In this paper, we study the curvature radiation in the 
rotating pulsar magnetosphere. The curvature emission of 
single particle and a bunch of particles streaming along the 
magnetic field lines are calculated using numerical simula- 
tions. The paper is organized as following. In Section 2, we 
present the details of our calculation for the emission from a 
single relativistic particle. In Section 3 the rotation distorted 
patterns of curvature emission are calculated for the emis- 
sion cone at a given height. The polarization profiles from a 
given height are calculated for three typical density models 
in the form of core, cone and patches for various emission 
geometries. The polarized profiles from the whole magneto- 
sphere are presented in Section 4. Our conclusions are given 
Section 5. 



2 CURVATURE RADIATION OF A SINGLE 
PARTICLE IN ROTATING PULSAR 
MAGNETOSPHERE 

The pulsar magnetosphere consists of relativistic plasma 
streaming along magnetic field lines. The magnetic field is 
generally thought to be an inclined dipole, which is given by 



B = B4 — f[3r{r ■ m) - m], 



(1) 



where B* is the surface magnetic field, _R* is the neutron star 
radius, r is the unit vector along r, and rh represents the 
unit vector for the magnetic dipole moment. Curvature radi- 
ation is generated when relativistic charged particles (elec- 
trons or positrons) streaming along the curved field lines, 
the electric field from a single particle is 

n X [{n — v) X a] 



E(t) = 



(1 



(2) 



with n the wave vector unity, v the particle velocity and 
a the acceleration. In previous studies of curvature radia- 
tion in pu lsar magnetospheres (JBenford fc Buschaueilll977l : 
iGangadha ra 2010), the rotation of magnetosphere is usually 
neglected, thus the particle velocity v is simply parallel to 
the tangential direction of the B field at the emission point, 
and the acceleration a is in the same B field line plane as the 
emission point. However, the co-rotation of a particle with 
the magnetosphere changes v and a, and affects the inten- 
sity and polarization of the particle emission significantly 
(jBlaskiewicz et al.lll99in . In this section, we discuss the cor- 
rection to the particle velocity and acceleration in detail, and 
give a precise equation to calculate the curvature radiation 
of a single particle in rotating pulsar magnetosphere. 



2.1 The particle velocity, acceleration and 
emission location 

In the co-rotating frame, the relativistic particles or bunches 
are streaming along the open magnetic field lines with a 
velocity of 
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Figure 1. The modified magnetic field direction, velocity, 
Lorentz factor and acceleration of the particles along the 
last open field line considering the co-rotation of pulsar 
magnetosphere. Fig. (a) shows the magnetic field direc- 
tions in the fixed lab frme (solid line) and the co-rotating 
frame (dashed line). In Fig. (b), (c) and (d), the solid 
lines show the precise relativistic velocity, Lorentz factor 
and acceleration given by equation (|4]) and ([7]), the dot- 
ted lines are approximations given by equations ([5]) and 
(l8t . Our calculation chooses the last open field line, which 
is in a plane writh azimuthal angle of —10°. The other pa- 
rameters involved in calculation are P = Is, a = 70°, and 
7 = 400. 



= v'b' 



(3) 



here b' — B/\B\ is a unit vector along the tangent direction 
of the magnetic field line, v' is the velocity magnitude in 
terms of light speed c. The particles also experience rotation 
with velocity Vi — fl x r/c. In the lab frame, the total 
velocity v should be the addition of v' with Vi by Lorentz 
transformation, which reads 

V = (4) 

'Jril + V.-V') ^ ' 

here v-r is the magnitude for the rotation velocity. Its cor- 
responding Lorentz factor is 7r = 1/Vl — v'r. Equation Q 
gives the precise particle velocity in the lab frame. The ve- 
locity V can also be written as 



Kb + Q, X r/c. 



(5) 



Here, k is the function of v and t)r, b denotes the 
magnetic field direction in the fixed lab frame, which 
is the Lorentz transformation of b' . The difference 
between b(solid line) and 6' (dashed line) is very 
small for a low emission height (r < 5007?*) as shown 
in F ig, [l] (a). Therefore, mo s t of the previous stud- 
ies (|Blaskie^vicz et al.l Il99ll : iKumar fc Gangadharal 
120120 can ignore the difference and use b ' inst ead of 
b in equation ((Sj. In iBlaskiew^icz et al.l (|l99lh . k is 
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Figure 2. The coordinates for pulsar emission geometry, ii rep- 
resents the rotation axis, m is the dipole moment. The magnetic 
inclination angle is a with respect to ii. Sight line has impact an- 
gle of 0, cutting the beam. The rotation phase is ifi = fit starting 
from the Cl — m plane. The magnetic polar and azimuthal angles 
are 9 and (h. 
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Figure 3. Tangential emission points (magnetic polar angle Bq 
and azimuth angle ipo) and the corresponding curvature radii p 
for various rotation phases, (p. The dashed lines represent the case 
without rotation, while the solid lines are the results by taking 
the rotation in to account. T he dotted lines show the results given 
bv lBlaskiewicz et al.l 11199 ll ). The parameters used for calculation 
are P = ls, a = 30°, /3 = 3°, 7 = 400, r = 250iJ*. 



limited to be a constant, since they just consider the 
emission from low heights where Vr is very small. To 
coordinate the total velocity with the light speed, 
iGangadharal (|2005l ) introduced the evolvement of k 
w^ith height by assuming the Lorentz factor to be 
a constant, as shown with the dotted line Fig. [TVc). 
However, according to equation (|4]), the Lorentz fac- 
tor for the particle is changing with height shown by 
the solid line in Fig.[l][c). The precise relativistic ve- 
locity given by equation (|4]) and the approximation 
are shown in Fig. [l](b). Within several hundred of 
Ri,, the difference is very small. 



The radiation of a relativistic particle is beamed in its 
velocity direction v. It requires the alignment between the 
wave vector unity n and v for the observer to receive con- 
siderable radiation. At a given emission height, the emission 
location {do, (j)o) at each rotation phase ifi (see Fig. [2]) can be 
determined by solving n = v, which is the tangential point 
of the trajectory along the line of sight. Here (^0, <^o) is the 
polar coordinate of magnetic axis , with (j)o = c orrespond- 
ing to the Ct-m plane. Blaskiewicz et al.l (119911 ) has done 
similar calculations by analytically solving n = Kb+flx r/c. 
Our numerical calculations develop their results further to 
higher magnetosphere using equation (|3)). The emission lo- 
cations from our calculatio n are shown in Fig. [3l w hich is 
very close to that given bv lBlaskiewicz et al.l (|l99l|). Rota- 
tion demonstrates its effect by shifting the emission locations 
towards an early phase. 

The particle acceleration is also affected by the co- 
rotation with pulsar magnetosphere. The charged particle 
streams along the magnetic field line, whose trajectory is 
described by r = re sin^ 6 in the co-rotating magnetic axis 
frame. Here r^ is the field line constant, 9 is the polar an- 
gle from the magnetic axis. The particle acceleration in the 
co-rotating frame reads 



V c ( 



(6) 



here b = dr/d9 is the magnetic field line tangent and its 
magnitude jbj — (re/\/2) sin 6*^/5 + 3cos(2S). The accelera- 
tion in the lab frame can be calculated by 



dv 
dT' 



(7) 



Substitute the approximation of v (equation^ in the above 
equation, we have 

AC c db hicdhi; „;■„,„ 

""7^9^ + 1619^'' + ^''"'''' + "''^"''''/''^' ^^^ 

w^hich is equivalen t to equation (7) in 
iKumar fc Gangadharal (|2012l ). Substitute the pre- 
cise velocity equation Q in equation Q, we can also 
obtain the analytic equation of a without approximation, 
which is too complicated to write done here. The numerical 
results for the precise acceleration from Lorentz transfor- 
mation and the approximated one given by equation ((8)1 
are shown in Fig. [Hd). The approximated treatment for 
the acceleration remains good only at low emission height 
of several hundred 7?*, but becomes significantly different 
from the precise acceleration at outer magnetosphere. 

The curvature radius of the particle trajectory in the 
lab frame is approximately given by p ~ a? ja. In the upper 
panel of Fig. [31 we calculated the curvature radii for the tra- 
jectories at each rotation phase using the acceleration. The 
leading side tends to have smaller curvature radius com- 
pared to the case without r otation, while the trailing side 
behaves in the opposite way. iThomas fc Gangadharal (|2007t ) 
pointed out that the rotation induces a significant curvature 
into the particle trajectories, which are deflected towards the 
rotation direction. The magnitude of the deflection on lead- 
ing side becomes larger compared with the no-rotating case, 
while the deflection on the trailing side is weakened. O ur re- 
sults is consistent with lThomas fc Gangadharal (|2007l ). The 
modifications of particle velocity, acceleration and emission 
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Figure 4. The circular path approximation (solid curve) to the 
actual particle trajectory (dashed curve). At the tangential emis- 
sion point, the particle velocity vg is parallel to the line of sight 
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Figure 5. The components of acceleration at each emission point 
parallel and perpendicular to the velocity direction. The solid 
lines are for the corrected acceleration including the rotation ef- 
fect, while the dashed lines are drawn without considering the 
rotation. The parameters used here are P = Is, a = 30°, /3 = 3°, 
7 = 400, r = 50R^,. 



location due to the co-rotation of particles with pulsar mag- 
netosphere have significant roles on pulsar radiation prop- 
erty, which will be discussed in Section 3. 



2.2 Curvature Emission from a single particle 

Knowing the particle velocity, acceleration and emission lo- 
cation of the relativistic particle, we can directly calculate 
its curvature radiation fields. Note that according to equa- 
tion ((5}, the emission location, where n • u = 1, is the point 
where the radiation electric field is the strongest. Before and 
after this point within the angle of I/7, the particle still 
emits in the direction of n but the magnitude of the elec- 
tric field becomes weakened (see Fig. |4]). The total electric 
field emitted from this particle should be the coherent inte- 
gration of the emission from the whole trajectory, especially 
near the tangential emission point. 

To study the emission from a relativistic particle, we set 



up a reference frame (eoi, 602, vo) at the tangential point 
along the ling of sight in the instantaneous trajectory plane, 
shown in Fig. U) Here the subscript "0" stands for the pa- 
rameters at the tangential point, vq and ao are the velocity 
and acceleration at the tangential point, vq is the unit vec- 
tor along Vo- The other two axes of the frame are described 
as 



€01 — eo2 X Vo, 
Vo X ao 



602 



Vo X ao 



(9) 



eoi is in the vo-ao plane and almost parallel to ao, and eo2 
is perpendicular to the plane. In Fig. (5] we compare the 
two orthogonal components of ao along ijo {a.o\\) and eoi 
(aox), and find |ao|||/|aox| ^ 10"^ for various particle tra- 
jectories at the corresponding pulsar rotation phases. This 
means that the contribution of ao|| t o the total em ission 
can be neglected. It was pointed out bv lJacksonI (119751 ) that 
the radiation emitted by an extremely relativistic particle 
subject to arbitrary accelerations is equivalent to that by 
a particle moving on an appropriate circular path. The ap- 
proximate circular path around the tangential point for the 
particle trajectory has a radius of curvature of 



Po = c- 



|ao| 



(10) 



here we suppose ao ~ ao_L, since |ao|| | <C |aox|. The velocity 
and acceleration of the bunch in the vicinity of the tangential 
point can be expressed by 

V = |do|[cos( )iJo + sm( jeoij, 

po pa 



ao|[-sm( )vo -f cos( JeoiJ. 

po Po 



(11) 



The position vector for the particle in the circular path reads 



r-ci, 



-pod, 



(12) 



with a — a/|a|. The radiating electric field of the relativistic 
particle at frequency u in direction n is the Fourier Trans- 
form of E{t), which is governed by 



E{u]) = 



qe- 



"Ro/c 



2ttcRo 



nx[{n-v)x a] ioj(t-n-r. 



■''"'dt,(13) 



here, Rq is the distance between the circular path center and 
the observer. 

We should note that the integration in equation 1131 
is highly osc i llatin g and not easy to integrate. Following 
iGangadharal (|2010f ). we expand the integrand to the third 
order around t = and then get the analytic expression of 
the integration. The total integration can then be calculated 
directly from the analytic equation. The detailed description 
of this method is given in the Appendix B of 'Gangadharaj 
(|2010l ). 

In our calculation, we use the assumption of a circular 
path to describe the trajectory of the relativistic particles 
around the tangential emission point. To demonstrate its 
feasibility, we further compare in Fig.[6]the integrand term in 
equation (|13[) near the tangential emission point using actual 
traject ory, similarly done by iKumar &i Gangadharal 
([2013) and circular path approximation. The magnitude of 
the integrand is highly oscillating when the particle is away 
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Figure 6. The integrand term in equation I I13I I near tlie tangen- 
tial emission point (wliere t = Os) considering the co-rotation of 
the particle with pulsar magnetosphere. Here A = |ti X [{n — v) X 
a]/{l—n-v)'^\, St = i— n-rcir/c, and Acos{io&t) are the magnitude 
of the integrand. The dotted lines are for the actual trajectory of 
the particle in the meridional plane, while the solid lines present 
the results using the circular path approximation. The left and 
right panels show the results for lower and higher emission radius 
(and corresponding frequencies), respectively. The other param- 
eters used for calculation here are P = Is, a = 10°, /9 = 3° and 
7 = 400. 



from the tangential emission point about |t| > 4 x 10~ s 
(or about lOi?*), so that the integration outside is almost 
vanished due to the oscillation. The central part near i = 
between the two dashed lines in the figure dominates the 
result of integration. In this central part, the difference by 
using the assumed circular path and the ideal trajectory is 
negligible. Therefore, the electrical fields calculated from the 
circular path approximation are accurate enough but much 
easier. 

The observed emission is the transverse component of 
the electric field perpendicular to the line of sight n. We have 
to transform the calculated E{uj) to the lab frame XY Z, 
with Z along the line of sight n, and fl in the XZ plane. 
The Stokes parameters are defined as follows 



/ = 


- Ex Ex + Ey Ey , 


-- 


ExEx^EyEy, 


u = 


= 2Re[E*xEy], 


V = 


= 2Im[ExEY]. 
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Figure 7. Pulse profiles of intensity and position angles emitted 
from a single particle at the tangential emission point at various 
rotation phase. The calculation is made for different heights r = 
lOOi?*, 200il* and SOOiJ*. The solid line presents the PA profile 
without rotation. The other parameters used here are P = Is, 
a = 30° and 7 = 400. 
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Figure 8. The fitted parameters Oflti /3flti ^'PO and 0pao vary 
with emission heights. The dots are the values by fitting the PA 
profiles shown in Fig. [7] with the rotating vector model (equa- 
tion [I5ll. The solid lines present the best fit to the dots, which is 
described by the equations given in each panel. 



2.3 Pulse profile of curvature emission from a 
single particle 

Previously, only the radiation exactly along the velocity di- 
rection of the relativistic particle is considered for the pulsar 
profile. It means that only the radiation by a particle at (6^0, 
(j)o) with a velocity towards the observer contributes to the 
observed intensity at a given rotation phase. Using equa- 
tion (I13fl . we can calculate the polarized emission from a 
single particle at various phases to get the intensity and po- 
larization profiles, as is sho wn in Fig. [T] The results are very 
similar to those given by iBlaskiewicz et al.l (|l99ll ). where 
they use the total emission of the particle in the emission 



cone at the tangential point instead of the integration along 
the trajectory. The emission is almost 100% linearly polar- 
ized since the emission at each phase comes from a single 
particle. The leading part of the intensity profile is brighter 
than the trailing side because of the rotation induced cur- 
vature for the particle trajectory. The most obvious are the 
polarization position angle curves shifted towards a later 
phase due to the aberration effect. 

The shift of a PA profile caused by rotation is larger at 
a larger emission height due to the relatively larger rotation 
velocity fl x r. To show the shifts varying with the emission 
heights, we can fit the shifted PA curves at different heights 
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by using the rotating vector model (|Koniesaroff et al.lll970l ) 



tan(0F 



sin(iy3 — A(/Po) sin a 



sm{a + P) cos a — cos(a + (3) sin a cos{ip — Aipo) 



(15) 



Here (^pao describes the PA profile shift along the per- 
pendicular direction in Fig. [71 Aipo represents the phase shift 
of a PA profile. Evolution of the four fitting parameters (ast, 
/3fit, Aipo and 0pao) with height are shown in Fig. (8] Ro- 
tation induced aberration not only delays the polarization 
sweep with positive Aipo, but also shifts t he entire PA swee p 
downward with negative (/>pao ((Hibschman fc Aronsll2001f ). 
Both of these shifts {4>paq and A(po) are the first order func- 
tions of r. Pulsar emission geometry is also distorted by the 
rotation, since the fitting values of aat and Pat change with 
heights. P A c urves arc steepened for /3 < but fiattened 
for /3 > (Isiaskiewicz ct al. 1991 ). The fitting value of the 
dipole inclination angle aat increases with r. The infiuences 
of rotation on am and Pat determination are both second 
order functions of r, as Pat — P <x r'^ and aat — a oc r^. 



3 CURVATURE RADIATION AT A GIVEN 
HEIGHT 

A relativistic particle streaming along a magnetic field line in 
magnetosphere will emit in a cone of I/7 around its velocity 
direction. Thus the observed emission at a pulsar rotation 
phase comes from the particles not only along one magnetic 
field line with n ■ b — 1, but also the nearby field lines 
within the I/7 emission cone. In this section, we will cal- 
culate the curvature radiation from the particles at a given 
height steaming along a bunch of magnetic field lines. 

3.1 Emission pattern at a given height 

The emission received at one rotation phase actually con- 
tains the contributions from all nearby field lines around the 
tangential emission point at a given height. Fig.[5]shows the 
emission patterns for the polarized emission with or without 
rotation at different heights and impact angle p. The pat- 
terns are plotted in the coordinates (6, <j)) around the mag- 
netic axis for different rotation phases. The central point 
of the patterns corresponds to the tangential point where 
the velocity aligns with the line of sight. The patterns ex- 
tend to approximately 0.15° in magnetic polar angle 9, and 
about 2°-4° in magnetic azimuth cj>, which agree with the 
1/7 emission cone since we take 7 = 400. The total inten- 
sity I distributes in an ellipse manner around the central 
maximum. The linear polarization L has almost the same 
patterns as I but with smaller magnitudes. The circular po- 
larization V is shown in two antisymmetrical lobes with -|- 
and — signs corresponding to the left and right hands, with 
the maximum intensity of V as about 20% of the peak total 
intensity. 

Our calculations in Fig. [9] make it clear that the rota- 
tion systematically distorts the emission patterns compared 
with the case without rotation in which the long axes of 
/ and L beams and the zero line of the V are all vertically 
aligned. Since particles at a larger emission height have com- 
paratively larger rotation velocity and hence trajectories are 
much more bent. For the emission from higher radii, the 



Table 1. The sign of the net circular polarization of particle 
emission for various density gradient at an emission height of 
50i?*. 



/3> /3< 



dN/ae > 

ON/ 89 < 



dN/d<f> > 
dN/d<t> < 



distortion of beam pattern due to rotation becomes more 
serious. The beam patterns for opposite impact angles (see 
Fig. [5^ and [S}i for ±/?) are almost the same except for the 
opposite signs in the V patterns. The beam pattern for a 
smaller impact angle P tends to have large distortion (see 
Fig. O), because the rotation induced acceleration is more 
dominating for the particles in the field lines near the mag- 
netic axis, compared to the curvature acceleration. Note that 
the beam patterns for the phases with opposite signs, i.e., 
ip — ±5° or ±10°, are asymmetric. The strongest distor- 
tion happens at a late phase {ip > 0°), because the rotation 
bends the trajectories towards the rotating direction. The 
long axes of the beam patterns for an extremely high emis- 
sion radii, i.e., r ~ 400-R* or higher, are rotated by more 
than 90° compared those without rotation. 



3.2 Polarization profiles for the emission from a 
given height 

The observed Stokes parameters at each rotation phase are 
the integration of all visible beam patterns from all particles 
in all field lines. It is generally believed that pulsar radia- 
tions are coherent, and the width of the coherent bunches 
in the 8-(j) plane should be much smaller than the radiating 
wavelength for the coherency to be efficient. The typical size 
of the 1/7 emission pattern at a given height is much larger 
than the coherent width. Therefore, the emission at a given 
height should be added incoherently, which reads 
eo+s9 ^4,0+64, 



Ir = / N{r, 

leo-se Jti>o-6</> 

Qr = I / N{r, 

Bo+se 1-4.0+64, 



I r sin 



e,(i)) Q r sinMM^, 



Ur 



Vr 



N{r;6,cj)) U r sin MM;/), 



'o+se /•4>o+S'l> 



N{r, e, 0) V r sin 



(16) 



Bo-SB J40-S4, 



Here N{r, 6, cj>) serves as particle density distribution func- 
tion, the half-width 56 and 5(j) of the emission pattern 
are chosen to be 0.1° and 2° to cover the pattern of the 
emission cone. For simplicity we assume that the den- 
sity distributions across r, 9 and (p are independent, i.e., 
N{r,9,(j}) = n{r)f{9)g{(j>). The density variation along the 
height n{r) oc r~^, and n{r) is a constant for a given emis- 
sion height. 

If the particle density, A^, keeps constant in all field lines 
at a given height, the observed circular polarization should 
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Figure 9. Emission patterns for the total intensity /, linear polarization L and circular polarization V around the tangential emission 
point calculated for various heights and /3 at different rotation phases (ip = 0°, ±5°, ±10°). The patterns without rotation are always 
similar for different r and /3, so that we just show one for r = 50-R* and /3 = 3°. The two axes for plots are the magnetic polar angle 8 
and azimuth angle <fi. The other parameters used for calculations are P = Is, a = 30°, 7 = 400 and u = Lu/2Tr = 400 MHz. 



be zero in all the cases since it is always symmetrical with 
the central point, and the integrations of particle emission 
from many field lines will smear out the circular polarization. 
However, in most cases the density distribution across field 
lines are not uniform. Modulation of TV across 9 and may 
result in net circular polarization. For example, in the case 
without rotation in Fig. ^jp, the pattern of V is symmetric 
with respect to </!>, so that only the density gradient along </!> 
can produce net V, while the density gradient along 6 can 
not. The detailed discussion for this case can be found in 
iGangadharal l|201(J ). In the cases of rotation, the observed 
beam patterns become complicated due to the distortion. 
The circular polarization depends on the density gradients 
along 9 and (j), following the rule shown in Table [l] for the 



emission from a typical emission height r = 50-R* (see Fig.|9^ 
and f). 



3.2.1 Cone-core model 

The observation of pulsar mean profiles suggested a classical 
pulsar beam model with a central co r e surrounded by one 
or multiple hollow cones (|Backeilll976l : iRankinll 19931 ) . These 
cone-core regions in pulsar magnetosphere are filled with 
relativistic particles streaming along the open magnetic field 
lines, whic h is generated by sparking proce ss on neutron star 
polar cap l|Ruderman fc Sutherland! 1 1975l l. We explore the 
emission profile for the particle density distributions in the 
form of core and cone in the polar cap region streaming 
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Figure 10. Polarization profiles (/, L, V and PA) at 400 MHz 
from a height of SOiJ* for the core density distribution. The solid, 
dash-dotted and dotted lines represent /, L and V , respectively. 
The density N , magnetic polar angle, 9, and the density gradient 
dN/dd vary vifith the rotation phase ip. The left and right panels 
are plotted for opposite impact angles /3 = ±1°. The other param- 
eters used here are P = Is, a = 30°, 7 = 400 and trg^ = 0.0005. 




Figure 11. Same as Fig. 1101 but for the density distribution in 
the form of a cone. Left and right panels stand for the emission 
profiles from 30 R* and 50 -R* , respectively. The other parameters 
used here are P = Is, a = 30°, ,9 = 5°, 7 = 400, ag^ = 0.0002 
and s = 0.7. 



out from the neutron star surface and extending to high 
magnetosphere. 

The core density distribution is written as 



f{e) = /oexp 

g{<P) = 1- 



2ai 



(17) 



Here /o is the maximura density, 6* is the polar angle of 
a field line footed on the neutron star surface in the polar 
cap region, ag^ is the characteristic width of the gaussian 
distribution on the star surface. The polarized profiles from 
particles distributed in this core at a height of r = 50i?* are 
shown in Fig. llUI Due to rotation, the "S" shaped PA curve 
is shifted towards a later phase, w hile the intensity curve t o 
an earlier phase, as predicted by iBlaskiewicz et al.l l|l991f ). 
Single sign of circular polarization is generated because of 
the density gradient along 6. In the case of positive /3 (see the 
left panels of Fig.[TOj, the left hand CP {V > 0) is generated 
because dN/d6 < for all phases, but the right hand CP 
for negative /3 (see the right panels of Fig. I10|l . Both of them 
agree well with Table [l] The core density model always has 
dN/d9 < and dN/dtj) = 0, so that the net CP always has 
a single sign for all rotation phases. 

The cone density distribution reads 



f{9) = /oexp 
ff(0) = 1- 



2ai 



(18) 



Here 0*iof = sin~''(-v/i?*/re_iof) is the polar angle maximum 
on the star surface corresponding to the last open field line, 
the dimensionless parameter s describes the peak position 
for the maximum density in terms of Otiof- The two sets 



of polarized profiles for the cone density model are shown 
in Fig. [11] with the same impact angle of /3 = 5° but two 
different emission heights. For a higher emission radius of 
r — 50i?* (left panels), the line of sight cuts through the 
density cone close to the magnetic axis, which results in 
double intensity components. For each component, dN/d6 
changes signs, which results in the sign reversal of circular 
polarization. The opposite dN/d6 variation causes the op- 
posite sign reversals of V. For a lower emission height of 
r — 30Rt (right panels), sight line will cut across the edge 
part of the density cone. Therefore only one component can 
be observed. The density gradient dN/dO < 0, so that the 
circular polarization has one sign across the pulse phases, 
very similar as the case of core emission. 

3.2.2 Patch model 

The patch density model describes the different density com- 
ponents randomly distributed within the open field line area. 
The patch foots located on the neutron star surface in the 
polar cap region can be described as 



f{e) = /oexp 



sO^, 



2ai 



9{<f>) 



ffoexp 



^'^l 



(19) 



Here go is the maximum density constant similar to /o , (f>p is 
the peak position of a density patch in the magnetic azimuth 
direction, and a^ is the corresponding characteristic width. 
We do not use the subscript "*" for <^ and (jip, because they 
do not change with the height. 

The position of a patch {s9i,iof, 4'p) affects the density 
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Figure 12. Sketch map of the discussed patches in the open field 
line region. The signs of "+" and "— " in the patches represent 
the signs of circular polarization at a height of 50 iJ* . See detailed 
profiles in Fig. 1131 and 1141 
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Figure 13. Same as Fig. 1101 but for the density in a patch at 
4>p = with r = 50R*. Profiles for six impact angles are shown. 
We plot dN/d(f)/{s'm9) here instead of dN/d<f>, although marked 
as dN/d<f> in the figure. The other parameters chosen here are 
P = ls, u = 400 MHz, a = 30°, 7 = 400, s = 0.7, ag^ = 0.0006 
and (T0 = 0.1. 



gradient along 9 and 0. For a given patch, different cuts of a 
line of sight will cause different profiles. Here, we calculate 
the profiles for typical patch positions in the open field line 
region, with various impact angles. Figure [12] is the sketch 
map of the positions of the discussed patches. 

Firstly we consider a patch located in the meridional 
plane with (j}p = 0. Sight lines cutting the patch with differ- 
ent /? give different polarization profiles, as shown in Fig. 1131 
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(e) <f)p = 225° 
Figure 14. Same as Fig |13l but for other patches in Fig. 1121 
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(g) </.p = 315° 
Figure 14. — continue 



For j3 — 2° , the dominating density gradient is ON/ 89, which 
is large and has a single sign at all the phases, and dN/dcj> 
is small. The circular polarization keeps to be of the right 
hand. When /3 increases to be 4°, there exits a little sign 
reversal of CP in the trailing side of the profile, mainly be- 
cause dN/dcj) is relatively strengthened. For the sight line 
cuts across the patch center of /? = 5.7°, ON/ 89 reaches 
its minimum, dN/dcj> dominates. The weak circular polar- 
ization shows a sign reversal. For large /3 > 5.7°, negative 
ON/ 89 dominates again gradually, which gives the reversal 
of weak CP and then the left hand CP. 

Now, we consider the discrete patches (see Fig. I12|l lo- 
cated at different positions around the magnetic axis with 
(/>p = 45°, 90°, 135°, 180°, 225°, 270° and 315°. The po- 
larized profiles from cutting through the inner, central and 
outer parts of the patch are plotted for each patch in Fig |14l 
Similar to the patch of (/>p = 0° , the V profiles always have a 
single sign, or the sign reversal, or the opposite sign for the 
line of sight cutting the inner, central and outer parts of a 
patch. For the patches at cj>p — 90° and 270° , we have V^ ~ 
when the line of sight cuts across the central part of the 
patch. The V profiles for the patch at (j>p = 180° has almost 
the same dependence on /3 as the patch of (pp = 0°, except 
for the reversed circular polarization. For other patches away 
from the meridional plane, the V profiles always depend on 
the variation of dN/d(j> and ON/ 86 as listed in Table [T] But 
for patches with cfip not close to 0° or 180°, the density gra- 
dient along (j) dominates the sign of circular polarization, 
because the rotation-induced distortions are weak. For con- 
venience, we mark the final sign of the V profile for different 
/3 in each patch of Fig. [121 

For a given patch, even when the impact angle /3 is 
fixed, the polarization profiles will be different for various 
emission heights. The polarized emission from two symmet- 
rical patches with a fixed /3 but different height (lOi?*, 20-R* 
and 40-R* , respectively) are shown in Fig. [151 Naturally pro- 
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Figure 15. Same as Fig |13l but for two symmetrical patches 
{4>p = ±30°) at heights of lORi,, 20/?* and 40_R*. The impact 
angle is fixed to be /3 = 3°. The other parameters used here are 
p = Is, u = 400 MHz, a = 30°, 7 = 400, s = 0.7, ag^ = 0.001 



and 



'<!>■ 



0.15. 



files gradually become wider with the increasing of height. 
They also gradually shift towards an earlier phase due to 
the aberration effect. The most important is the evolution 
of the circular polarization with emission heights. At a low 
height of lOi?*, two components have different signs of CP, 
and each with one sign only. At a height of 20i?* (the mid- 
dle panel), each component has the sign reversal of CP. At 
a height of 40-R* , each component has one single sign of CP 
again without reversal but the opposite hand compared to 
those at 10 i?*. This is mainly because at different heights, 
the sight lines with a fixed P cut through the inner, middle 
and outer parts of the patches, which should have different 
circular polarization due to the dominating of dN/dcj). 



4 EMISSION FROM OPEN FIELD LINE 
REGION 

In general, a dipole field is assumed for pulsar magne- 
tosphere. The tangential points of the field lines of 6 ~ 
2/3 6nm can produce emission detectable for a given line 
of sight n. Here 9nm is the angle of the sight line from the 
magnetic axis m, and 9nm = P when the sight line goes to 
the meridional plane. Curvature radii are different for the 
tangential points at different heights, so that characteristic 
emission frequencies are different. This is the basic princi- 
ple for the "radius to frequency mapping". However, the 
emission from one height actually has a wide-spread spec- 
trum, the characteristic frequency is only the frequency of 
the peak-intensity. 

In fact, one receives the emission from the whole radi- 
ation region as long as the line is of S ~ 2/3 9nm in the 
open field line region (see Fig. I16|) . The retardation effect 
should be considered for the emission from different heights. 
The emission profile from the whole open field line region is 
given by 



I{ip) = Ir{ip + r/ric)dT, 



Qiv) 



n(v3-|-r/ric)dr. 
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Figure 16. Emission regions for different density models. Tlie 
thick line outlines all the tangential emission points of field lines 
of 6 ~ 2/30Tim for a given line of sight, with 9nm as the angle 
between wave vector n and magnetic momentum m. The grey 
areas represent the emission regions for the core and cone/patch 
regions. 



t/(^) 



Ur{(p + r/nc)dr, 



V('fi) = / K-(<^ + r/nc)dr 



(20) 



Here Ir, Qr, Ur and Vr are given by equation (|16p . and r/ric 
is the retarded phase for the emission from different heights. 

The emission from a patch in the meridional plane with 
s = 0.7 is plotted in Fig. I17t i. with the polarization gener- 
ated at different phases and heights. The peak intensity Ip 
(given by equation [16} varies with the heights. The signifi- 
cant contribution to the integrated emission profiles in the 
bottom panel comes from the region of 8-20i?* . For r < 8i?+ 
and r > 20i?* , the emission from this patch points away from 
the observer (see the sketch map of Fig. fTS)) . The dominating 
contribution comes from the peak height around r ~ 12i?» 
for the maximum of plasma density (see the dashed line in 
the left panel of Fig. ll7b .V For the patch with a higher and 
wider emission region of s = 0.3 fFig. 117b ). the dominating 
emission comes from the height of about 130-R*. There the 
aberration and retardation effects are rather strong, and the 
phase for the peak intensity and CP are shifted to an early 
phase. 

We plot the emission profiles from the peak height in 
Fig. 1181 which are very similar to the integrated total inten- 
sity and polarization profiles. The emission from the peak 
height always has a sign reversal of the circular polarization, 
because the density gradient along the magnetic azimuth (p 
has the sign reversal. The aberration and retardation effects 
are very weak for r — 12i?* and the depolarization effect is 
also negligible. 

We noticed that the polarization profiles from a patch 
are almost the same for different frequencies. If the patch is 
wide and the density gradient is small, e.g. the almost uni- 
form distribution, the net circular polarization from a patch 
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Figure 17. Total intensity and circular polarized emission of the 
curvature emission from two patches, plotted for various phases 
and heights. The linear polarization is very similar to the total 
intensity / except for a slightly smaller magnitude. The integrated 
emission profiles are plotted under each map. The variations of 
the peak intensity Ip (solid line) and Nr^ (dash-dotted line) are 
also plotted against heights in the left. The pulsar parameters 
used here are the same as those in Fig. 1131 




(a) r 




(b)r 



Figure 18. The emission profiles from the peak heights 
two patches in Fig. 1171 



should be very small. Note however that the peak height is 
different for different frequencies, and the integrated circular 
polarization may change for different frequencies. 

For the core density model, the emission region of 
6 — 2/39nm does not have an upper limit (see Fig. I16p . 
The intensities from various emission heights, up to thou- 
sands of i?*, are still significant compared with that from 
the peak height of r ~ 2507?* (see Fig. 119b ). It is related to 
the integrand term Nr'^ in equation (|16p . which determines 
the final value of Ip. The density variation can be written as 
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Figure 19. The peak intensity varies with the heights (solid lines) 
for core density model (panel a), and for cone model (panel b). 
The density variations, Nr^, are also plotted with the dashed 
lines. The parameters are the same as those in Fig. 1101 for core 
and Fig. Ill I for cone. 
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Figure 20. Same as Fig. 1171 except for the curvature emission 
from core (panel a) and cone (panel b). Lorentz factor damping 
along the particle trajectory is assumed as eQuatior l22l with an 
initial Lorentz factor at the star surface of 70 = 1000 and the 
damping factor K-y = 0.01. The other parameters used for cal- 
culation here are the same as those in Fig. 1191 



Nr^ = n{r)f{e)g{(j))r^ oc 



exp 



R. 



2ai 



(21) 



and plotted by the dashed line in Fig. 119b .. The total emis- 
sion profile is very viride in phase because of the retarded 
emission from large heights, and the polarization will also 
be vanished. 

In a practical case, the damping of the Lorentz factor 
of a particle along r is natural. For example, an exponential 
decay of the Lorentz factor was introduced by iQiao et al.l 
l|2001l ) to analyze the inverse Compton scattering process in 
the magnetosphere. Here we take the slow variation of the 
Lorentz factor with height as 



7 = 70 exp(-KT. 



R. 



R* 



(22) 



where 70 is the energy for the secondary particles near 
the neutron star polar cap, and k^ is the damping fac- 
tor, and we calculate the total emission profiles as shown in 
Fig. 120b .. The emission region is constrained to be between 
20i?* and 1207?*. The total polarized profiles resemble to 
those from a height of 507?* (Fig |10|) . The circular polariza- 
tion contributed from each height has the same sign, so the 
final total circular polarization has only one single sign. The 
total intensity and polarization profiles are asymmetry due 
to the retardation effect. 

For the cone density model, the density Nr'^ and the 
peak intensity vary with height as shown in Fig. 119b . The 
emission can be produced even at a very large height. The 
dominating emission comes from about 30 R* to 300 7?* . The 
profiles become very wide (see Fig |20b ) because of the strong 
retardation effect. With the damping of the Lorentz factor 
of particles as the core model, the total emission profiles and 
7p(r) are shown in Fig. 120b . The cone density only produces 
one profile component, instead of two. Because the emission 
from low heights are very strong, and two components from a 
higher regions are weak and only add the outer wings in both 
sides of the profile from the lower height. The total circular 
polarization is very weak also, because of the depolarization 
of the emission from different heights. 



5 DISCUSSIONS AND CONCLUSIONS 

In this paper, we depict the detailed scenery for curvature 
radiation in rotating pulsar magnetosphere for the first time. 
The rotation affects the velocity and acceleration of particles 
and hence the emission property. Circular path approxima- 
tion for the particle trajectory is used at the emission point 
to calculate the radiating electric field (see equation I13p . 
The emission from particles at the tangential points of field 
lines and the nearby field lines in the I/7 emission cone is 
considered for a given phase and height. The polarization 
profiles from a given height and the whole open field line 
region are calculated for three possible density distributions 
in the form of core, cone and patches. We find the following 
conclusions: 

(i) Rotation not only shifts the PA curves along the rota- 
tion phase, Aipo, but also causes an offset of the curve, 0pao, 
both of which are the first-order functions of the emission 
height. Its influences on a and /? determination follows the 
second order functions (Fig. [S}. 

(ii) Rotation distorts the patterns for the I/7 cone more 
seriously if it originates from larger height and/or smaller 
impact angle (Fig. [9)l . The density gradients across the pat- 
terns will result in the net circular polarization. 

(iii) For the patch density model, the V proflles from the 
same height can have a single sign or the sign reversal or the 
opposite sign, depending on where the sight line is cutting 
(Fig. ll3lT^ . The polarization profiles resulting from a fixed 
impact angle also vary with emission heights (Fig. I15|l . 

(iv) The central peaked core will have the circular polar- 
ization of only one hand, which is usually significant and 
always has a single sign (Fig. 120b). The cone emission at 
a large height is separated into two components with the 
opposite sign reversals for circular polarization. But at a 
lower height there is only one component with a single sign 
of V. The total circular polarization is very weak because 
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of the depolarization of the emission from different heights 

(Fig.EOb). 

Our calculations for the curvature radiation from par- 
ticles distributed in the form of core, cone and patches 
can be compared with the observed pulsar polarization 
profiles. The sign reversal of circular polarization are al- 
ways observed in the centr a l component of the mean profile 
jLvne fc ManchesteJll988l : iRankinI 1 19831 : ICil fc Snakowskil 
1 199(3 ). iHan et alj (19981) pointed that there also exists sign 
reversal of circular polarization for the cone components for 
some pulsars. Such observational facts can be explained if 
the emission comes from central or non-central patches. The 
CP with a single sign may come from the core emission. 

Here we did not consider the possible propagation ef- 
fects, i.e., the wave mode coupling effect ( Wang ct al. 201CI ) 
in the magnetosphere. Our calculations can explain some 
pulsar polarized emission, but it can not match various ob- 
servation facts. It should be noted that we consider just the 
curvature emission from particles with a single 7 and as- 
sume the simple distributions of particle density. The ac- 
tual energy and density distributions of particles in the 
magnetosphere are not clear at p resent. There e xist other 
emission mechanisms , like ICS (|Xu et al.l I200C1 ). plasma 
process (JBeskin et al.lll988l '). Moreover, we treated the co- 
herency in a simple way by assuming the coherent bunch 
to be a huge point charge. However, the actual coherent 
manner, bunc h size and shape affect the radiation pattern 
IjBenford fc Bu schaucr 1977). We used the vacuum dipole 
magnetic field as the basic model for our calculation, but 
there exists the possible distortion of the dipole magnetic 
field d ue to rotation, which has not been taken into ac- 
count (jDyks fc Har ding 2004). The magnetic field structure, 
such as t he globally twiste d, self-similar, force-free magne- 
tosphere (|Thompson et al.ll20 02l. can also lead to the mod- 
ification of the radiation pattern. All these will make the 
actual emission pattern much more complicated than the 
simple calculations presented in this paper. 
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